The Covariance Function of a Simple 

Trunk Group, with Applications 

to Traffic Measurement* 

By V. E. BENES 

(Manuscript received July G, 1960) 

Erlang's classical model for telephone traffic in a loss system is con- 
sidered: N trunks, calls arriving in a Poisson process and negative expo- 
nential holding times; calls which cannot be served at once are dismissed 
without retrials. Let N(t) be the number of trunks in use at t. An explicit 
formula for the covariance R(-) of N(-) in terms of the characteristic 
values of the transition matrix of the Markov process N(-) is obtained. 
Also, R( ■ ) is expressed, purely in terms of constants and the ''recovery" 
function, i.e. the transition probability Pr\N(t) = N \N(0) = N};R(-) 
is accurately approximated by R(0)e ni , with n the largest negative char- 
acteristic value, itself well approximated (underestimated) by — E{N(-)}/ 
R(0). Exact and approximate formulas for sampling error in traffic meas- 
urement are deduced from these results. 

I. INTRODUCTION 

A theoretical study of sampling fluctuations in telephone traffic meas- 
urements is useful both in designing procedures for measuring traffic 
loads and in interpreting field observations. Hayward 1 and Palm have 
given an approximate formula for the sampling error incurred when 
observations of the numbers of calls in existence are made at fixed in- 
tervals of time. Their formula has the disadvantage that it is derived 
for a probabilistic model (of the traffic) in which there is an infinite 
number of available trunks. Thus there is no limit to the number of 
calls which can be in progress at one time, and no congestion. Two im- 
portant parameters, the number N of trunks in the group, and the prob- 
ability p N of loss, are left out of account. For this reason the practical 

* This work was completed in part while the author was visiting lecturer at 
Dartmouth College, Fall-Winter, 1959-GO. 
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application of this model is usually restricted to large groups of trunks 
which are lightly loaded. 

In this paper we derive and study the covariance junction of the sim- 
plest stochastic model of a finite group of N trunks. The sampling error 
in traffic measurements can be calculated exactly from the covariance. 
We find formulas for the magnitudes of fluctuations of observed traffic 
for both periodic and continuous observation. The exact formulas lead 
to simple approximations similar to Hayward's, which take account of 
the number of trunks. Our results are summarized and discussed in 
Section II. 

We shall use A. K. Erlang's classical probabilistic model for a group 
of trunks, described as follows: 

i. Holding times of trunks are mutually independent, each with a 
negative exponential distribution. Time is measured in units of mean 
holding time. 

ii. Epochs at which calls arrive form a Poisson process of intensity 
a > 0, independently of the holding times. The offered load is then a 
erlangs. 

iii. There are N < cc trunks; calls which find all .A/ of these trunks 
busy are "lost," and are cleared from the system. 

These assumptions determine a Markov stochastic process N(t), 
— co < t < co , the number of trunks in use at time t. N(-) is a random 
step-function fluctuating in unit steps between and N. As is well known, 
N(-) has stationary probabilities \p„ , n = 0,1,- • •, N] given by the 
(first) Erlang distribution 

a" 
n\ 

Vn = ~N~—k 

T°L (1) 

= equilibrium probability that n trunks are busy. 

With this choice of absolute probabilities, N(-) is a strictly stationary 
process, whose mean and variance are respectively 

mi = a(l — p N ), 

a = mi.— ap N (N — Wi). 

The probability p N of loss is shown in Fig. 1, the fractional occupancy 
A r_1 mi in Fig. 2, and the variance a 2 in Fig. 3. 
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Fig. 1 — Probability pn of loss. 
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Fig. 3 — Equilibrium variance. 



II. DISCUSSION, SUMMARY AND CONCLUSIONS 

The covariance R(t,s) between samples N(t),N(s) of the stochastic 
process N(-) is the average of the product of N(t) and N(s), minus 
the product of the averages: 

R(t,s) = E{N(t)N(s)} - E{N(t))E{N(s)}. 

Since N(-) is a stationary real process, we have R(t,s) = R( \t — s\). 
The function R(-) is called the covariance function of the process iV( • ). 
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It can be written as 

B(t) = lim E{N(t + u)N(u)\ - E[N(t + u)\E{N(u)} 

" N /A V ^ 

= E wp„.E n¥i{N(t) - n|tf(0) = m} -(E mp*.) , 

„»=0 n=0 \ro=0 / 

where {p m } are the stationary (or equilibrium) probabilities given by 
(1), and 

Pr{AT(f) - n | iV(0) = m) 

denotes the transition probability that n trunks are busy at time t if m 
were busy at time 0. The function R( • ) expresses the average depend- 
ence or correlation between samples of JV( • ) taken at times t apart. 

The principal practical use of the covariance function R( • ) in the 
theory of telephone traffic is in computing theoretical estimates of sam- 
pling error incurred in traffic load measurements. Two methods of meas- 
uring traffic, the switch-count and the time-average, are considered in 
this paper. In the switch-count, n observations [xi , • ■ • , x n , Xj = N(jr), 
j = i f . . ., n \ of the random process are made at intervals r apart; the 
average 

-I>Or) mi£ x ,-n*8. 

n j=i n ;=i 

is then used as an estimate of the carried load mi = a — ap N . This 
method is important economically because it is cheaper to scan trunk 
groups periodically than to observe them continuously. The number t 
is the scan interval, and the number S„ = Xx + • • • + x n is called the 
(total) number of paths in service, in n observations. Table I lists actual 
mean holding times, scan intervals used and resulting values of t for 

Table I — Holding Times, Scan Intervals and Values of t 



Type of Call 


Typical 

Holding Time 

(seconds) 


Scan Interval 
(seconds) 


Ratio t of Scan 

Interval to Holding 

Time 




U.S.A. 


Europe 


U.S.A. 


Europe 


Local Call 
Long Distance Call 
Originating Register Holding 
Time 

No. 5 Marker Holding Time 


100-200 

200-600 

10-15 

0.25-1.0 


100 

100 

10 or 100 

10 


36 
36 
36 


1 to i 
Jto* 
1 to 1 

or 
10 to 7 


J to J 

g to is 
4 to 2 
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various types of call. The variance of n~ 1 S n is expressible in terms of 
the co variance R( • ) as 

Var{,fX} = n* % (n - \j\)R(jr). (3a) 

In the time-average, the continuously recorded sample average 

M(T) = T~ l f N(t)dt, 

•'o 
is used to estimate the carried load. The variance of this estimate is 

Var{M(T)} = 2T~ 2 I (T - t)R(t) dt. (3b) 

Jo 

Thus the mean square error of both these methods of traffic measure- 
ment can be calculated theoretically if the co variance R( • ) is known. 

In formula (2) the co variance function is expressed in terms of the 
stationary probabilities \p n \ given by the Erlang distribution, and the 
transition probabilities 

p mn (t) = Pr[N(t) = n\N{0) = m\. 

In the theory of telephone traffic, the particular transition probability 

pss{t) = Vv\N(t) = N | AT(0) = N\ 

has been singled out (in Refs. 3 and 4, for example) as a suitable "re- 
covery" or "relaxation" function that is characteristic of the dynamic 
behavior of the Markov process N(-) in point of the undesirable "all 
trunks busy" condition. 

We shall show that a much more cogent reason than this can be ad- 
duced to support the importance of the recovery function to traffic 
theory: the covariance function R( • ) can be expressed entirely in terms 
of the recovery function and the offered load a. In other words, a single 
one of the (JV + l) 2 transition probabilities appearing in formula (2) 
suffices for determining the covariance function, and this one is the 
recovery function j) N n{ • ). This fact is a theoretical justification of the 
intuitive view that the recovery function is important, for now the vari- 
ances of w _1 aS„ and of M(T) are expressible using only the recovery 
function. 

We next give a summary of the contents of the remaining sections; 
this is followed by an account of specific results and conclusions. 

An exact formula for the covariance R( ■ ) is stated and discussed in 
Section III, and derived in Section VII. The formula readily yields a 
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rigorous upper bound which appears to give a close approximation to 
R(-) itself. In Section IV the recovery function p NN ( ■ ) is given, and it 
is shown how the covariance may be expressed in terms of the recovery 
function by a convolution integral. The variance of n -1 <S„ is studied in 
Section V; an exact formula, and an approximating upper bound [based 
on the upper bound for R( ■ )], are both obtained. The variance of the 
time-average M(T) is considered in Section VI; again, an exact formula 
and an approximating upper bound are found. 

The covariance function R(-) is bounded from above and closely 
approximated by a single exponential function 

R(t) ^ *V 1 ', a' = E(0), n < 0. 

Here 

tr 2 = R(Q) 

= equilibrium variance of N( • ) 

= (load carried) — (load lost) (average number of idle trunks), 

and the reciprocal time constant )\ in the exponent is the dominant* 
characteristic value of the "rate" or "transition" matrix of the differ- 
ential equations satisfied by the transition probabilities. Alternately, 
ri is the root of least magnitude of a Poisson-Charlier polynomial. The 
root n is shown as a function of offered traffic a for N = 1, • ■ • , 8 in Fig. 
4, and is tabulated in Table II. 

A lower bound for r\ , depending only on the mean and variance of N(-), 
is derived in Section VIII by making use of the fact that the matrix of 
the differential equations for the transition probabilities is symmetriza- 
ble. For low values of offered traffic per trunk, i.e., a/N < 1, this bound 
can be used to approximate >i . In any case, the bound is a convenient 
starting place for the use of Newton's method. The bound is the ratio 
— nh/cr", which satisfies the inequality 

7 ^ n < - 1, 

with 

mi = equilibrium mean of N( • ) 

= load carried = a(l — ps), 

a 1 = equilibrium variance of N(-) 

= (load carried) — (load lost) (average number of idle trunks). 

The approximation n ^ — m^o is illustrated in Fig. 5. 

* I.e., that of least magnitude (among the nonzero characteristic values). 
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I 



-V 




Fig. 4 — Negative of the root n of smallest magnitude as a function of load a 
for # = !,-••, 8. 



Table II — Negative of Dominant Characteristic Value ri 


a 


N = 4 


N = 5 


.V - 6 


N - 7 


N = 8 


1 


1.043967 


1.011448 


1.002421 


1.000421 


1.000062 


2 


1.249464 


1.112166 


1.045044 


1.015806 


1.004800 


3 


1.582363 


1.326321 


1.172257 


1.084025 


1.037229 


4 


2.000000 


1.629624 


1.383389 


1.222707 


1.121762 


5 


2.477548 


2.000000 


1.663799 


1.427870 


1.265214 


6 


3.000000 


2.422137 


2.000000 


1.689991 


1.463798 


7 


3.557618 


2.885474 


2.381627 


2.000000 


1.710891 


8 


4.143703 


3.382497 


2.800900 


2.350437 


2.000000 


9 


4.753426 


3.907677 


3.251918 


2.735363 


2.325514 


10 


5.383178 


4.456828 


3.730121 


3.150052 


2.682770 
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Fig. 5 — Illustration of the approximation n = — mi/a 2 . 



By the "infinite trunk" model we shall henceforth mean the stochas- 
tic model for telephone traffic determined by all the same assumptions 
that we made in the Introduction, except that N = °o ; i.e., an unlimited 
number of trunks is postulated. Riordan 5 and BeneS 6,7 have considered 
this model; Hayward 1 based his sampling error formula on it. 

It is widely believed that the "infinite trunk" model is applicable to 
large groups of lightly loaded trunks. Such a belief is gratuitous until 
comparisons with a model having a finite number of trunks are made. 
Studying the covariance function of the simple finite trunk group en- 
ables us to make some of the needed comparisons; e.g., the variances of 
n^S,, and M(T) in the two models are of particular interest. Knowledge 
of the covariance R(-), however, is also relevant to the other three cases 
to which engineers are loath to apply the "infinite trunk" model, viz.: 

i. large groups of heavily loaded trunks, 
ii. small groups of lightly loaded trunks, 
iii. small groups of heavily loaded trunks. 
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The variance of n~ S„ is bounded from above and approximated by 
the formula 

( l — ~ 2nX 1 

VarfrT 1 ^,} ^ n~V <ctnh X - —— csch 2 X , (4) 

where n is the number of observations, and 

X = — — = — \ (scan interval) (dominant characteristic value). 

The exact formula for the variance of n~ l S„ in the "infinite trunk" model 
is 

n~ a sctnh \t — — = csch 2 \ r > . (5) 

The upper bound (4) for the finite group is compared with the exact 
formula (5) for the "infinite trunk" model in Fig. 6, which shows each 
formula as a function of the scan interval t for various n, for a = 20 
erlangs offered to 20 trunks. The curves suggest that the upper bound 
for Var|w~ S n \ for N < <» is consistently less than the corresponding 
variance in the "infinite trunk" model. As might be expected, increasing 
the scan interval r improves accuracy for the same number of observa- 
tions. This is because the covariance function is positive, and monotone 
in | / |. 

The variance of M(T) is bounded from above and approximated by 

Var{M(7 1 )} ^ 2<x 2 *- - ~^~ nT , (6) 

where T is the length of the time-interval of continuous observation, 
and a 2 and r x are, as before, the variance of N( • ) and the dominant 
characteristic value, respectively. The exact formula for the variance 
of M(T) in the "infinite trunk" model is 

2atL^±l. (7) 

Since 7\ < —1, and a is always less than a if N < °o, the "infinite 
trunk" model overestimates the variance of M{T) if applied to a finite 
group. This conclusion is illustrated in Fig. 7, which shows the formulas 
(6) and (7) for a load of 20 erlangs offered to 20 trunks. For an observa- 
tion time of 10 mean holding times the "infinite trunk" formula (7) 
applied here would overestimate the variance by about 500 per cent. 
This is about as extreme a case as would occur in practice. Fig. 7 also 
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Fig. — Comparison of variance of S„/n for finite and infinite trunk models. 



depicts a "mixed" formula obtained by replacing a by a 2 in the "infinite 
trunk" formula (6); for 10 mean holding times the "mixed" formula 
only overestimates the variance by about 100 per cent. Thus most of 
the discrepancy is due to the difference between a and a. 
Our conclusions are set down in the following list: 

1. The average dynamic behavior of the process N(-), as described 
by the covariance function R(-), can be adequately determined from 
the dominant characteristic value n and the variance a 1 of N( •)• 

2. The same parameters, n and a, suffice to give simple approximat- 
ing upper bounds for the sampling error incurred in both periodic and 
continuous observation of N( • ). These bounds depend on the size N of 
the trunk group. 
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Fig. 7 — Comparison of variance of M(T) for finite and infinite trunk models. 



3. In terms of r% and a it is possible to check the applicability, for 
theoretical estimates of sampling error, of the "infinite trunk" model 
which assumes N = °o . 

4. The "infinite trunk" model, applied to finite trunk groups, con- 
sistently and often grossly overestimates the sampling error. The over- 
estimation occurs largely because o- 2 is always less, and for heavy traffic 
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is much less, than a, the (Poisson) variance of N( • ) in the "infinite 
trunk" model. 

5. In terms of n and a it is possible to design sampling procedures for 
traffic measurement that depend explicitly on the number N of trunks 
in the group. By these methods, a given accuracy can be obtained with 
less observation, and thus at lower cost, than the "infinite trunk" model 
would require. 

6. Hence for finite groups of trunks traffic sampling procedures which 
are based on the "infinite trunk" model tend to be wasteful, particularly 
for heavy traffic. The parameters r t and <r 2 provide a systematic way of 
tailoring the measurement procedure to the number of trunks in the 
group. 

III. THE COVARIANCE FUNCTION 

To state the formula for R( ■ ) we need the "sigma" functions* de- 
fined (see Riordan ) as 

f \ a 

with m (but not /c) a nonnegative integer. These functions are connected 
with the Poisson-Charlier polynomials 

Pn ( X ) = «" /8 (»0*E/-i) n "{")i!^. 

by the relation 

a k (m) = (-a*) m (ml)"*p«(-*). 

(SeeRef. 10, p. 33.) 

For fixed N and a, let ri,r 2 , •■ • ,r N he (in order of increasing magni- 
tude) the N zeros in the variable s of the polynomial tr,+i(N). In Sec- 
tion VII the covariance is shown to be given by (the exact formula) 

Bit) = -a 2 V s t ,f* 5 II [1 - (rj - r.-r 1 ] (8) 

where p N is the probability of loss. It has been shown f that the zeros 
Tj are all real, negative, and distinct; all are less than —1, and consecu- 

* The a notation is copied from unpublished work of H. Nyquist . The functions 
themselves were introduced into traffic theory by Palm. 8 

t The earliest reference appear to be Haantjes 11 in 1938. See also Ledermann 
and Renter. 12 
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tive pairs are separated by at least unity. Fig. 8 shows these roots for 
N = 1, 2, 3 as functions of a. 

Now rj is always negative, and the terms of the product satisfy 



1 - 



1 



>0; 



(9) 



Tj - Ti 

hence the sum in (8) has all terms negative, so that 

R(l) ^ 0, all t. 

The correlation between successive samples is thus always positive. It 
is obvious from (8) that 
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Fig. 8 — Roots of the first three a-i unctions. 
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Since r x is the root closest to zero, the value of (8) is only increased 
if the Tj in the exponents of (8) are replaced by n . Using (9) and (10), 
we conclude that 

^ erV 1 ' - fi(0 = $(<), 00 

where 

««) = « 2 p, e r . n , _ )t n ' , 

and 

9-1 , c-> / ■ . . \— 2 ^. 2 47r" — 30 _f 

f(0 ^ a-p w e £ ( 3 + 1) ^ a Vs ki e 

£ ( 0.3933 )ap N e~ l . 



The approximation R(t) = crV 1 ' is illustrated in Figs. 9 and 10. It ap- 
pears to be fairly accurate, especially for light loads. 

The upper bound oV 1 ' for R(t) should be compared with the rigorous 
formula (see Riordan 5 and Benes" ) 

R(t) = ae~ , 

which holds for the "infinite trunk" model. In this model the equilibrium 
distribution of occupancy is Poisson, so that 

R(0) = a 1 = Varf7V(0l = E{N(t)\ = a, 

and the "time constant" of the exponential is unity, since time is meas- 
ured in units of mean holding time. 

The difference between the "infinite trunk" model and the "finite 
trunk" model in point of the covariance can be understood by con- 
sidering the effect of congestion, which is present in the latter. Conges- 
tion affects the upper bound formula most directly through the value 
of the variance a 2 . It is obvious intuitively, and borne out in Fig. 3, 
that as a increases a 2 must eventually decrease to zero. This behavior is 
not mimicked by the "infinite trunk" model, for which a 2 = a. 

The finitude of N, i.e., congestion, affects the bound <rV r| (| in two ways: 
(a) the "time constant" is not unity but the smaller number — (fi)" 1 , 
so that the rate at which dependence between samples of N( • ) de- 
creases (as a function of the interval between samples) is larger than 
in the "infinite trunk" model; this "time constant" decreases as the 
traffic a increases, because, as illustrated by Fig. 4, n is a monotone 
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Fig. 9 — The covariance R(t) for N = 5 trunks, a = 10 erlangs, with the ap- 
proximate formula It(t) ~ aV 1 '. 

decreasing function of a; (b) the value of R(0) ( =o- 2 ) is not a but the 
generally much smaller number 



2 

a = 



= a(l- P ,)[l-ap.( a -^--l)] ) 

= a[l - p N (l + N - a + ap N )]. 

The last form shows that a 2 < a for all a and N. In fact, it is obvious 
intuitively that 

<t 2 = Wi — ap N (N — mi) < nil < a. 

A simple approximation for the dominant root n can sometimes be 
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used to make the approximation R{t) = o-V ri<l more useful. It is shown 
in Section VIII that 



/Ml 



carried load 
load variance 



^nj 



i.e., —mil a is a rigorous lower bound to n . Fig. 5 suggests this bound 
gives a fairly good approximation to r\ if a/N < 1. Hence a simple 
approximate formula for R( • ), valid for a/N < 1, is given by 



,, . N f carried load . 

= (load variance) exp< — = — = : t 

load variance 



(12) 



We know that R(t) ^ oV 1 " 1 ' 1 and that —mi/ a- 2 < n ; hence replacing 
fi by —mi/ a' tends to correct the error in the upper bound formula. 
The formula (12) is illustrated in Fig. 11. 
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Fig. 11 — Comparison of R(t) with(r 2 e" (m «'» u for N = 8 trunks, a = 4 erlangs. 
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IV. THE COVARIANCE IN TERMS OF THE RECOVERY FUNCTION 

It has been shown 3 that the Laplace transform of p N s( • ) is 

Vs(N) 

By expansion in partial fractions we find that 

Pnr(i) = Ps ~ E^- n(l " ZT^fJ' 
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(13) 



The sum assumes only negative values, and so p NN ( ■ ) decreases mono- 
tonically to the loss probability p N . The recovery function is illustrated 
in Fig. 12. 
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Fig. 12 — Recovery function for N = 5 trunks, a = 10 erlangs. 
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We now observe that for each j — 1, • • • , N, 

f (t - u) e -<'-^» du = f-=*l _ J±L . (i 4) 

^ (rj + l) 2 r j + 1 

By comparison of formulas (8) and (13), and use of (14), one finds that 
R(t) = a*p N [ (t - u)e- u - u) [p NN (u) - Ps\du + «rV + Cte~ l , (15) 
where 

c=-a> P „± i n(i--^-). 

i-irArj + 1) ,vy \ r,- - n/ 

This formula expresses #( • ) in terms of p/f N ( ■ ) by a simple convolution. 
To evaluate C explicitly we note that 

C = -r 2 - T °"m-i(N — 1) a -i 



2 I 0" 

ap 4a 



+ s)<r 8+1 (N) 1 + s_ 



where a_i is the first coefficient in the power series expansion of the left- 
hand term in the bracket. One finds 

_ <ro(N - 1) N 
a - x <ro(N) a ' 

C = ap N (^ - 1 + P") 

= ap N (N - mi) 

= (load lost) (average number of idle trunks). 

V. THE VARIANCE OF THE NUMBER OF PATHS IN SERVICE 

We assume that n observations [xj ,j — 1,- • •, n] of JV(-) are made 
during an interval of equilibrium, so that 

Cov{x it Xj) = R(\ i - j | r), 
where r is the scan interval. Then with 

S„ = X\ + X2 + • • • + x„ 

= number of paths found in service, 
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we find that 

Var{£ n ) =i'{EE x t x\ - E 2 JE ft] 



= EECov^vr,) (16) 

= E (n-\j\)R(jr). 

j——n 

To give an exact formula for Var{ S„] we note that 
f; e - 2lmlu = ctnh u, 

mm— oo 

E | m I <f ,,m| " = -# E e- 2 " ,u = | csch 2 u 

m— — oo ttM m=l 

and 

00 oo 

E/ \ —tutu 1 -2nu T~* I ,v„ I „- 2 l m l" 

(m — n)e = \e £^ \ m \ e 

m=n m=— oo 

= i<T 2 " u csch 2 w. 
Then also 

£ ( n -|j|) e - 2|/,tt =71 E «" , W- E Me~ 81 

+ 2 E (m-n)e- 2mu (17) 

/ , — 2nu \ 

, , (1 — C ) ,2 * 

= n ctnh M ^ csch M." 

Since the covariance B( • ) is a symmetric function given by (8), it can 
be seen that 



Var{n _1 5 B | = 






--v»S ' v -% (1 +V v ' /J (18) 



* Use of this identity was suggested by unpublished work of J. W. Tukey to 
which the author had access. 
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This formula is exact, given the assumptions. It is easily shown from 
formula (17) that the exact formula for the variance of n _1 »S n in the 
"infinite trunk" model is 



i *a <ctnh |t — — csch 2 %t> , 



illustrated in Fig. 6. 13 

Returning to the case of finitely many trunks, we can obtain ap- 
proximating upper bounds to formula (18) for Var{n _1 ^„} by using the 
results of Section III on the covariance function. It can be seen from 
the arguments leading to (17) that replacing the roots r,- by ri in the 
hyperbolic functions in (18) increases the values of the expressions in 
square brackets; this replacement is equivalent to using the upper 
bound 



2 r\t 

a e 



for R(t) in formula (16). Hence 



-2nX 



Var{rr\S„} ^ n"V jctnh X - L_i — cscn 2 ,\ . ; i<n 



where 



X = — — - = — ^ (scan interval) (dominant characteristic value). 

Since a e ri is close to R(t), we may expect that the overestimate (19) 
gives a good approximation to the actual variance. This approximation 
is conveniently plotted as a function of X for various n in Fig. 13. 



VI. THE VARIANCE OF TIME AVERAGES 

It follows from formulas (3b) and (8) that 

Vary N(t) dt\ = c + Cl T + o(e~ T ) (20) 

as T — * * , where 



„n(i--M 



«.-^P„g rHl + rtr 
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Fig. 13 — Upper bound to VarjS.I/wa*. 

is a negative constant, and 

Note that c and c x differ only in the power of r> that occurs in the 
denominators. The third term of (20) is positive; is given by 

equals — c at T — 0; and is of smaller order than e~ T because n < — 1. 
To evaluate c x explicitly, we note that 



C\ 



_ 2 p g.+i(^ - 1) _ l-p* _ a-i _ «-2 I / 9 , 
^° PAr U(l + «)Wtf) « 1 + a (1 + «) 2 J.-o' v- 
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where a_ 2 , a_i are respectively the first and second coefficients in the 
power series expansion of the leftmost term in the bracket of (21) ; these 
are given by 

<r (JV - 1) N 
a_ 2 = 



ao(iV) a ' 

1- -i 



dx xa x+ i(N) 



= N i 1 ~P\ 
a ap N 

To find Ci we must compute 

|" <r s+1 (N - 1) _ 1 - p N l = _d_ a x+1 (N -1) 1 
,.o L*(l + s)V 8+1 (iV) s _ " dx (1 + xy<r x+l (N) Jx=o ' 

This equals 



or 



(1 - Vn ) h - ^ log , x+1 (tf - 1) + A bg cr I+1 (A0 j 
Now the generating function of the <r-functions is 

*(s,z) = EgV.(n) = (1 - z)~V 2 
so that 

as n w — 1 1 

«fe Jx=o j-o (w — J)<T\{n) 

It follows that 

d = 2a 2 py[a_ 2 + a_j + (1 - Pa0(2 - $ S -i + £*)] 
= 2c7 2 + 2a~p N (l - p N ){\ - ^_x + fr) + 2aATpjr. 
The constant c can be evaluated in a similar fashion. 
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From the bounds (11) for R( • ) we conclude that 
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0^2. 2 ( e -^i^)-Var{M(T)} 

^ (0.3933)a*p* ( ' - ~ T \ + T ) . 



and since R(t) = aV 1 ', we may expect that the overestimate 



2 6 



n T 



- 1 - nT 



TW 



(22] 



is a good approximation to the variance of M(T). This approximation 
has the same form as the exact formula (7) for the "infinite trunk" 
model, because in both cases a single exponential is used for R( ■ ) in 
formula (8). The overestimate (22) is depicted graphically in Fig. 14. 
It was convenient to plot the ratio 

VartMCT)} 

<7 2 

as a function of the single parameter 

M = Ti T = —(dominant characteristic value) (observation time). 
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A simpler form of (22), valid for a/N < 1, results when we replace r\ 
by its lower bound 

mi carried load . 

«r load variance 

This replacement decreases the value obtained, i.e., moves the approx- 
imation in the direction of Var{M(T)\. 



VII. DERIVATION OF THE COVARIANCE 

The transition probabilities 

p mn (t) = Pr{N(t) = n | N(0) = m\ 
of N( • ) satisfy the Kolmogorov equations 
p mn (0) = 8 mn , 

j. PmN = ap m(A r_i) — Np mN , 

at 

d (23) 

Tl Pmn = (w + l)p m( „ + D +ap m( „_i) - (a + n)p mn , < n < N, 

at 

d 

-r PmO = Pml — apmo . 

Multiplying the nth equation by n, and summing on the index n, we 
find 

j t E{N(t) | N(0) =m\ = -E[N{t) \ N(0) = m} + a[l - p mN (0\, 

whence 

#{#(<) | N(0) = m) = me~' + a [ e~ (l - u) [l - p mN (u)] du, 

Jo 

By formula (2), the co variance is then 
B(t) = E mp m E{N(t) \N(0) = m\ - m* 

nt N 

= m i e~ t + ami(l — e~ l ) — m? — a I e~ (t ~ u) JZ mp m p mN (u) du, 

Jo m-0 
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where 



N 

Mi = 2 n% P» 



71-0 

for i = 1, 2 and {p n \ are the stationary probabilities given by (1). In 
particular, 

Wl = a(l - V s), (24) 

a = (m 2 - mi 2 ) 5 = [mi - apAr(iV - mi)]*. (25) 

The Laplace transform of 

Pr{N(-) = iV|JV(0) = m} 

has been determined 3 to be 

a N ~ m m\a.(m) 
msa. +l (N) " 

Therefore that of R( • ) is 

(s) = e R{t) dt = — — + -rp 

Jn 1 + S S(l 



/? 



+ «'«(! + «) 



a N ~ m mla.(m) 



(26) 



2 wp„ 



By ( 1 ), the last term of (26) is 

- ,d+^.w S w " (m) - 

It has been shown 9 that the "sigma" functions satisfy the recurrences 
<r s (m) = a,+i(m) - a.+i(m - 1), (27) 

m<r t (?n) — aa s (m — 1) + sa 8+ i(m — 1), (28) 

so that 

S N-\ N-l 

53 ma,{m) = a^2 <r*(k) + sXl v*+i(k) 
m— I *=o *.=o 

= aa s+1 (N - 1) + s<r, +2 (JV - 1), 

and 

<r, +2 (AT - 1) = N aa a+ i(N - 1) 
a. Tl (iV) s+ 1 (s + l)<r.+iW" 
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The foregoing identities yield the following simplified formula for R*(s) : 
R*( ) = m2 + ami aNp N 



1 +S 8(1 + «) (1 + s) 2 

_ ap N I" (t, + i(N — 1) "[ _ m? 
From (27) we find that the partial fraction expansion 



tr. 



AN - 1) _ v _ *r,(N)N\ 



Z 



<r a+1 (N) pi ( s - rj) II (rj - r t ) 

-£(«->■,)- n ,f 7 * 7 * ■ 

i=l ir*j Tj — Ti 

is valid, where [ryj are the zeros of <t s+ i(N). 

By a similar argument, since pn = <to(N)/<ti(N), 



Cg+l 



(AT - 1) 1 - p N N 



s{\ + 8)<r.+i(N) « «(1 + *) 



-V 



il 1 TT r i — 1 - r < 



7-i r, 1 + rj &j rj - r< 

Hence formula (29) can be inverted to give, for / ^ 0, 

R(t) = mi.e~ + omi[l — e~ ] — aNp N te~ — wii 



, /•' _ {( _ w) r jv _„ a /*■ 

Mi r y - TV J 



(30) 



2 _j . 2 »- -< 2 r> e T ' TT r i — 1 — U 

... +a P „Kc -,p,g- n - j _ ) - 1 g —- , 

where 

<r 2 = ra 2 — Wi 2 = equilibrium variance, 

and 

K = V I 1 TT r i ~ 1 ~ r » 
«ri(l + ri)»tt r y -r, ' 
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To evaluate K explicitly we observe that 

r (r. +l (iv - D_ _ _a_i_ fl^_i (31) 

where a_2 , a_i are respectively the first and second coefficients in the 
power series expansion of the leftmost term in the bracket of (31). Thus 

K = a_i + a_ 2 — 1 + ps ■ 
Now 

(l +S )V 8+1 (iV)- (,+S) .oW + u + s) 

fd a I+ i(iV- 1) 1 ^ (a - fj)" 1 tj rj - 1 - r, 

'|_<fa ^ +1 (.V) _U-i "*" pi (1 + ry) 8 M r ; - n 

From the recurrence (28) for the a-i unctions we find that 

differentiating with respect to .s and setting s = 0, we obtain 

d a, +1 (N - 1) 1 _ 1 ( _ <n(N - 1) \ _ _ 1 - Py 

a_1 ds it.+i(N) U-i a \ <*>(#) / ap N 

Clearly, 

«r„(JV - 1) iV 



0_2 = 



a (A7) ■ a 

and so 

1 — p N . N . . 
apAr a 

a"p N K = —a'. 

Thus the formula (30) for the covariance function R(- ) simplifies to 

R(, » = -4^Hw- (32) 

VIII. APPROXIMATION TO THE DOMINANT CHARACTERISTIC VALUE 

The differential equations (23) can be written in the form 

j Pit) = QP(D, 

at 
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where P(t) is the matrix of transition probabilities [pmn(t)} and Q is 
the matrix of the "transition rates": 



Q = 



-a I 

a (-a - 1) 2 

a (-a - 2) 3 





















; 


• 


N - 1 







a (- a - N + 1) 


JV~ 




a 


-N 




, • • • , tn • We define 







N „} 



M „ = I =n\a- n Z^, n = 0,l,--,N, 

Pn ;=oj! 

and we introduce an inner product for the space L 2 (ju) of (N + 1)- 
tuples of complex numbers by the definition 



(x,U) - 2 XnpnMn ■ 



N 

— > 



The matrix Q represents a symmetric operator on Li(fi), i.e., 

(Qx,y) = (x,Qy), x,y£Li(n). 
It is easily seen that 

N , N 

E--Zp. -i. 

,i=0 fi n "=0 

Qp = 0, for p = (p , pi , • • • . p N ), 
Qnfii = Qafij, i,3 = 0,\,---,N. 
The last identity implies that 

(Qx,y) = -sZ (yon - ViVj) — (m**< - «*i)i 

2 ,\; m.- 

(Qs,z) ^ 0, 

so (as we already know) all characteristic values of Q are nonpositive, 
being of the form (Qx,x) for some .t£L 2 (m)- 



(33) 

(34) 
(35) 
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From the extremal properties of the characteristic values of symmetric 
operators (e.g., Zaanen, 14 p. 383, Theorem 3) we conclude that 

n = max(Qx,x), 

the maximum being over all z£L 2 (m) which are not identically zero, 
and satisfy (x,x) = 1, (x,p) = 0, p being the vector of stationary 
probabilities, as in (34). 

We can now estimate r\ from below by choosing an appropriate 
vector x. We choose 

x n = - 1 , n = 0,\,---,N, 

crfi„ 

where mi and a are the mean and standard deviation of N{ • ) in equi- 
librium, given by formulas (24) and (25) respectively. Clearly, (x,x) = 
1 and (x,p) = 0, and 

if— l / I i „ \2 

(r . v v^ in - mi n + 1 - m,\ 

{Qx,x) = -a 2^ Pn [ ) 

n=o \ <t a / 

a(l - p N ) 



_mi < , 
a 

(See Kramer. 15 ) 
This approximation is illustrated in Fig. 5. 
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